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Summary 


A study of the effects of four wing-leading-edge modifications on a general 
aviation aircraft's stability and control parameters is presented. Flight data from 
the basic aircraft configuration and configurations with wing modifications are 
analyzed to determine each wing geometry's stability and control parameters. The 
parameter estimates and aerodynamic model forms are obtained using the stepwise 
regression and maximum likelihood techniques. The resulting parameter estimates and 
aerodynamic models are verified using vortex-lattice theory and by analysis of each 
model's ability to predict aircraft behavior. Comparisons of the stability and 
control derivative estimates from the basic wing and the four leading-edge modifica- 
tions are accomplished so that the effects of each modification on aircraft stability 
and control derivatives can be determined. 
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INTRODUCTION 


Unexpected entry of general aviation aircraft Into stall and spin flight condi- 
tions and the unsatisfactory handling characteristics that result in these regimes 
are major factors in aircraft accidents and fatalities (ref. 1). As a result of the 
work done through the NASA General Aviation Stall/Spin program (ref. 2,3), it was 
indicated that modifications of an aircraft wing's leading-edge promise improved 
stall/spin characteristics. The previous stall/spin research program included 
results from spin-tunnel tests, free-flight radio control model tests, full-scale 
aircraft flight tests, and rotary balance wind tunnel tests. Though various wing 
modifications were shown to improve the stall/spin characteristics of general avia- 
tion aircraft, the majority of this work has been qualitative in nature. In general, 
these studies either determined modes of instability associated with each type of 
wing modification, or have determined some of the basic aerodynamic characteristics 
of each geometry. A complete flight data analysis of the effect of wing-leading-edge 
modifications on an aircraft's stability and control parameters is of assistance in 
quantifying which of the various wing-leading-edge modifications would be most favor- 
able for improving the stall/spin charcteristics of such aircraft. Where the previ- 
ous work was less systematic in nature, this study is an effort to systematically 
develop complete aerodynamic models for these wing configurations. 

In an effort to analyze the effects of wing-leading-edge modifications on an 
aircraft's stability and control derivatives, a study was completed for a full-scale 
general aviation aircraft on which wing-leading-edge modifications were placed. 
Flights were flown in the pre- and post-stall regimes with the angle of attack in the 
range of between approximately 4 and 24 degrees. The pre-stall regime is defined for 
this paper as the angle of attack range where 1) the lift curve is nonlinear and 2) 
a is less than the static stall angle of attack. Post-stall will refer to flight 
conditions in which a is greater than the static stall angle of attack. The data 
that resulted from the flight tests was analyzed using the stepwise regression and 
maximum likelihood techniques. From these analyses, a complete set of stability and 
control derivatives were estimated. Though there was no wind tunnel data available 
for comparison, the parameter estimates were verified by using vortex-lattice theory 
to estimate theoretical stability and control values, and by inspection of the air- 
craft behavior prediction capabilities of the estimates. The stability and control 
derivative estimates from each wing geometry were compared and the effects of wing 
modification on the stall/spin characteristics were determined. 

Presented in this paper is a description of the flight test program that was 
performed, and of the data analysis techniques that were used. Also discussed are 
the aerodynamic models which were used in the parameter estimation techniques. A 
discussion of the stability and control derivative estimates is presented along with 
concluding statements and remarks. 


FLIGHT TEST AND DATA COLLECTION 

As a continuing project of the NASA General Aviation Stall/Spin program, a 
variety of general aviation aircraft possessing wing-leading-edge modifications have 
been flight tested in an effort to improve the stall/spin characteristics of such 


aircraft. Data from the flight tests of these aircraft have been collected and made 
available for analysis by NASA personnel and other researchers working in the area of 
stall/spin improvement. The data on which this study is based were provided by the 
staff of the Flight Dynamics Branch of the Low-Speed Aerodynamics Division at the 
NASA-Langley Research Center. 

The flight test program was carried out for a low-wing, single-engine general 
aviation aircraft. A three view drawing of the test airplane with the basic wing is 
shown in figure 2. The four different wing-leading-edge modifications were created 
by adding a drooped section to the basic wing. The resulting modif ications are shown 
in figure 3. They include a wing on which the full leading-edge span was modified by 
addition of the drooped section (FLE modification), a wing in which a gap was placed 
in the drooped section, creating a segmented leading-edge (SLE modification), and a 
wing in which the outboard sections of the span were modified by the droop and a 
fairing which smoothed the leading-edge from the drooped section to the basic wing 
(FOLE modification). A wing geometry similar to the FOLE modification but without 
the faired section was also analyzed. This irregular leading-edge geometry is known 
as the outboard leading-edge (OLE) wing modification. The physical characteristics 
of all the aircraft test configurations are presented in table 1. 

The test aircraft was fully instrumented to record angular rates, linear accel- 
erations, control deflections, angles of attack and sideslip, true airspeed, engine 
RPM, attitude angles, and altitude. The linear accelerations were measured by three 
accelometers mounted orthogonally near the aircraft CG location. Angular rates were 
measured by three separate rate gyros which were themselves orthogonally mounted. 

The rudder, aileron, and horizontal stabilator control deflections were measured by 
position transducers placed near those surfaces. Angles of attack and sideslip and 
true airspeed were obtained from flow direction and velocity (FDV) sensors, which 
were mounted on tip booms in front of each wing's leading-edge. The engine RPM was 
measured via a magnetic pick-up mounted on the engine flywheel. Attitude angles were 
measured by a three-axis miniature attitude reference system (MARS). The roll and 
pitch angles are measured with respect to a vertical position of the spin axis. 
Altitude was sensed by a pressure transducer in the aircraft's static system. 
Measurements were collected from the instrumentation and digitized at a rate of 
20 samples/sec and then stored on magnetic tape. 

The flight tests were executed such that the collected data covered a large 
range of attack and sideslip angles. Longitudinal, quasi-steady, acceleration- 
deceleration maneuvers were flown along with longitudinal, lateral, and combined 
large and small amplitude transient maneuvers. The quasi-steady maneuvers were 
maneuvers in which there was almost no short period excitation (q * 0) as the pilot 
very slowly pulls back the stick until a maximum controllable a is achieved, then 
slowly returns the stick forward until the original trim a is reached. From such 
maneuvers, one can derive lift and pitching moment curves. 

Combined maneuvers involve a mixture of both lateral and longitudinal excitation 
in the same test maneuver. The longitudinal motion was excited from trimmed condi- 
tions by the pilot applying simple stabilator doublets and pulses, while the lateral 
motion was excited from trimmed conditions by pilot applied sequential aileron-rudder 
or rudder-aileron inputs. For all the configurations tested, the aircraft was flown 
with a full aft CG location (31% c). 

The digitized flight data collected from each test configuration was corrected 
to reflect values referenced to the aircraft CG location and to reflect previously 
determined instrumentation calibration. 
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The flow direction and velocity (FDV) sensor measurements were corrected for 
wing upwash and sidewash. These corrections were determined from flight and wind 
tunnel tests. The velocity measured from the FDV sensors was corrected to reflect 
calibration runs flown with the test aircraft. The airplane was flown at various 
speeds over a measured ground trace. The true velocity as determined by the time 
required to cover the ground trace was compared to that measured from the FDV sensors 
and an appropriate linear calibration was determined. This calibration was then 
applied to obtain corrected airspeed. The angle-of- attack measurements were 
corrected by analysis of steady-state flight test data. The corrections were 
extrapolated to higher attack angles through the use of unpublished data obtained in 
the NASA-Langley 12-foot tunnel. 

Many of the measurements from the aircraft instrumentation were not made at or 
near the aircraft CG location. Therefore, angles of attack and sideslip, velocity, 
and linear accelerations were subsequently transformed to the aircraft CG location. 
The left and right boom measurements from the FDV sensors were transformed to the 
aircraft CG and then averaged together to obtain corrected airspeed and angles of 
attack and sideslip. The offset of the accelerometer package from the CG was known, 
hence, the accelerometer measurements were transformed to the CG location. The 
attitude angle measurements from the MARS gyro package were transformed to Euler 
angle measurements based on the angles at which the gyros were uncaged. The trans- 
formation and calibration equations used in the data correction are given in 
appendix A. 

In an effort to confirm that there were no bias errors in the boom or accelerom- 
eter measurements or scale factor errors in the boom measurements, the extended 
Kalman filter-smoother, fixed-point smoother of reference 4 was applied to the data. 
If the filter-smoother routine indicated measurement bias/scale factor errors, the 
data were then corrected so that compatibility between the various measurements 
existed. 


PARAMETER EXTRACTION TECHNIQUES 

The successful parametric modeling of the aerodynamics for an airplane operating 
near stall angles of attack consists of two phases. First, the aerodynamic model 
structure must be determined and second, the associated aerodynamic parameters (sta- 
bility and control derivatives) must be estimated. Though the aerodynamic model 
structure is known to be linear at low angles of attack, the existence of nonlinear- 
ities at higher angles of attack has been made known in several recent reports on 
flight test results, (refs. 8,17,20). 

Various methods exist (refs. 5,6,7,22,23) which allow for the determination of 
aircraft stability and control derivatives from flight data. Two of these methods, 
stepwise regression (SR) and maximum likelihood (ML), were used in this study. Each 
method differs in its level of simplicity, optimization criteria, and ability to 
determine aerodynamic model structure. 


Stepwise Regression 

Since a large number of possible nonlinear terms could contribute to the aero- 
dynamic function, some method must be developed that examines only the useful terms 
while ignoring those that are superfluous. One possibility is to look at all com- 
binations of linear and nonlinear terras. However, the number of models to be 
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considered grows too fast with the number of possible terms for such a technique to 
be practical. The use of the stepwise regression was suggested in reference 25. 

The stepwise regression examines each term as to its usefulness in improving the 
model (by reducing residual variance). 

Stepwise regression (ref. 8) belongs to a family of estimation techniques known 
as equation error methods. Equation error methods, which are based on the least 
square principle, estimate parameters such that the sum of squares in the difference 
between the measured aerodynamic forces and moments and those calculated by the 
estimated parameters is minimized. 

In the development of least squares methodology, assume that a dependent 
variable, y(t), can be closely approximated as a linear combination of independent 
variables. For an aircraft, the aerodynamic model can be formulated as 


y(t) = e 0 + 6^(0 + 


+ 6 


n-1 


x n-l (t) 


( 1 ) 


The dependent variable y(t) represents the aerodynamic force and moment coefficients, 
0} to 0 n -i represent the stability and control derivatives, 0 q is the value 
of the aerodynamic force or moment coefficient when the independent variables are 
zero, and x^(t) represent the n-1 independent variables, which are composed of the 
state and control variables or their combinations. 

At any particular time, t^; i = 1,2,...,N, the measured data can be used to 
approximate equation (1) by 

y(i) = 0 q + QjXjd) + . . . + e n _ 1 x n _ 1 ( i ) + e (i) (2) 


Because N observations were made, N linear equations approximating equation (1) 
result. These approximations contain an equation error term, e(i). Given that the 
number of observations is greater than the number of unknowns, N > n, the vector of 
unknown parameters, 0 , may be estimated using the method of least squares (ref. 5): 

0 = (X T X) -1 X T Y (3) 

where 0 is the n x 1 vector of parameter estimates, X is the N x n concatenated 
matrix of measured independent variables, and Y is the N x 1 vector of measured 
dependent variables. 

The properties of the least squares estimates depend upon the assumptions 
postulated about the equation error and the independent variables. Given the 
assumptions : 

1. e is stationary with zero mean 

2. e is uncorrelated with X 

3. X is deterministic 
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o 

4. e(i) are independent with mean zero and constant variance 

5. e(i) is normally distributed 

then the covariance matrix of parameter estimates will have the form 

E{ (0 - e )(8 - 0) T } = o 2 (X T X) _1 (4) 


where the variance of the equation error, 



can be estimated by 


s 


2 


as given by 


2 

s 


1 

N-n 


I e 2 (i) 
1 


(5) 


e(i ) = y(i) - y(i) 


( 6 ) 


given that the dependent variable can be estimated as 


yd) - e Q + e^d) + • . • + VjVi ( i) 


(7) 


It should also be noted that in real situations, the previously mentioned 
assumptions are not, in general, met. Thus, the least squares estimates can be 
asymptotically biased and can have higher standard errors than those predicted by 
equation (4). 

Before the least square concept is used to formulate the SR method, it will be 
of importance to define the random variable F which has an F-distribution with 
degrees of freedom n-1 and N-n as 


F 


m 

8 X Y - Ny 
(n-l)o 2 


( 8 ) 


o 

Another important variable is the squared multiple correlation coefficient, R , which 
is found by the relationship 


Am m O 

„2 . 8 1 X i Y - Ny z 
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Y Y - Ny 


(9) 
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where the mean value of y Is found by 


y -■£ I y(i) 

1 


( 10 ) 


The squared multiple correlation coefficient can be related to the variable F by the 
equation 


F 



(ID 


Also of consequence is the test of the significance of individual terras in the 
regression (partial F-test). A particular parameter, 0j, is deemed significant if 
its partial F-test criterion, Fp, given by 


F 

P 



* 2( V 


( 12 ) 


is greater in value than the tabulated F-distribution, F(nj, h2» a p ) » which possesses 
degrees of freedom hj * 1, - N-n, and a significance level otp (ref. 24). 

The stepwise regression procedure is a least squares estimator in which indepen- 
dent variables are inserted into the model, one at a time, until satisfactory results 
are obtained. The independent variables are inserted into the regression model based 
upon the partial correlation coefficient of the variables not currently included in 
the regression equation. Starting with a postulated regression model as in equa- 
tion (2), the first independent variable for the postulated model is chosen as the 
one which is most closely correlated with y. The correlation coefficient, r^y, is 
given by the equation 


■jy 


» * 

Jy 


1/2 


(13) 


which is composed of the sum of square terras defined as 


s jy = ^ “ x j ]( y (i) ~ y> 


( 14a) 
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(14b) 
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(14c) 


s yy " l ty(i) " y]2 


and where the mean value of is defined as 


x . 
J 


= £ l *<<*> 


(I4d) 


The independent variable, x^ , with the maximum r^ value is first selected for use 
in the model. If Xj is selected as Xp then the model is fit to the data by the 
expression 

y " e o + e i x i O 5 ) 


A new independent variable, Z 2 , is created by finding the residuals of X 2 after 
regressing it on Xp This is done by defining 

z 2 = Aq + Aj Xj (16) 


Then Z 2 is given by 


z 2 = x 2 - Aq - A x (17) 


The variables Z3, z^,..., z n _^ are found in a similar way b ^ regressing the vari- 
ables x 3 , x^,..., x n-1 on xj. A new dependent variable, y , is represented by the 
residuals of y regressed on Xj . Using the model given by equation (16), t.e.. 


* 

y 



■ Vi 


(18) 


a new set of partial correlations which include the variables y , Zj, ..., z n _i 
are formulated. Such partial correlations of z^ and y are related to the model 
containing the variable Xj . The next independent variable added to the regression 
model is one of the remaining independent variables x^ whose partial correlation is 
greatest. Again, the same process is completed and the correlations between the 
residuals of x^ regressed on the inserted independent variables and the residuals 
of y regressed on the inserted variables are found. 

When the partial F-test criterion of equation (12) is applied to the stepwise 
regression technique, model structure determination can be completed (ref. 8). As 
regression models are formulated, each variable inserted into the regression model is 
tested using the partial F-test criterion. If the variable does not meet the 
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criterion, that variable is dropped from the model. Hence, only statistically 
significant terms exist in the model, and the structure of the model is revealed. 
However, it must be noted that the partial F-test criterion can include too many 
terms in the model (overparameterization). Therefore, other statistical quantities 
should be taken into account when selecting adequate models. Some of these criteria 
are: 

1. maximizing F 

2 

2. maximizing the R coefficient 

3. maximizing F p terras 

4. creating a residual sequence which should appear uncorrelated and 
Gaussian 

Because all these criteria cannot simultaneously be met in practice, it then becomes 
a case of judgement and experience in determining which model is most adequate. 

Because the least square estimates are not explicitly functions of time, the 
independent and dependent variable measurements can be arranged in any order. This 
concept can be used to partition the measurements. Each partition is created by 
collecting the measurements into groups based on the values of one or more of the 
independent variables. Regressional techniques can then be applied to each parti- 
tioned set separately. The resulting parameter estimates for each partition then 
reflect the relationship of the dependent variable to the independent variables for 
that particular partition. Such a method is useful when large quantities of data are 
available, and if a statistically useful number of observations are contained in each 
partition. 


Maximum Likelihood 

Another method which is well suited for parameter estimation is the maximum 
likelihood method. In its simplest form, the maximum likelihood (ML) technique is an 
output error estimation method. Using such methods, parameters are estimated such 
that the sum of squares in the differences in the measured and calculated outputs is 
minimized. Whereas the SR method estimates the unknown parameters for each separate 
equation describing the dependent variable, the ML method must indirectly estimate 
the unknown parameters by solving a multivariate system. 

In the formulation of this method, it is assumed that only the measured outputs 
are corrupted by noise. The state, output, and measurement equations are, 
respectively 


X 

II 

£ 

/'-’N 

rt 

+ BU(t ) 

(19) 

Y = CX(t) 

+ DU( t) 

(20) 

Z i = Y i + 

V ± , i=l ,2 , . . . ,N 

(21) 
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where X, U, Y, and Z are the vectors of the aircraft states, inputs, outputs, and 
measured outputs respectively, A and B are the system and control-distribution 
matrices for the state equation, and C and D are the state-distribution and control- 
distribution matrices for the output equation. It is assumed that the measurement 
noise has the properties 

E{V ± } = 0 ; EtV^} = R6 

where R is the measurement noise covariance matrix and is initially an unknown. 

The maximum likelihood cost function (ref. 6) for the system can be written as 

J(0 ) = j l (Z t - y 1 ) T R' 1 (Z 1 - Y i ) +| In |R _1 | (22) 


where the notation, A , will always refer to the current estimate of the associated 
variable. R can be estimated by minimizing the cost function with respect to R. 
This minimization produces 


R - 


- y V V T 
N 4 Vi 


(23) 


Equation (23) is substituted into equation (22) and then the cost function can be 
written as 


J(0) = I(Z i - Y 1 )T*-l(Z l - Y^ ) + constant 


(24) 


Minimization of equation (24) with respect to the unknown parameters represents 
a nonlinear estimation problem. There exists several methods for its solution. The 
most common method used in aeronautical applications is based on the modified Newton- 
Raphson algorithm (ref. 9). In this algorithm, it is assumed that the output can be 
approximated as 


A 

A A Ay 

Y(6 ) = Y(0 O ) + A0 

0Q 


(25) 


where 0g is some initial estimate of the unknown, and A0 is the difference 
between the actual value of the unknown and the initial estimate of the unknown. 
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The cost function is minimized by differentiation with respect to 9, which 
yields 


3 J (0 ) 
30 


0 = -l gJr l v + l G R *G 40 

i i i 1 1 


( 26 ) 


where 



(27) 


are the sensitivity functions and the residuals, v, are given by 


v i = \ - ^(Oq) 


(28) 


Equation (26) is solved for A0 by writing it in the form 

a® - [X g[r"V 


(29) 


The sensitivity function G must be found by integrating the sensitivity equa- 
tions. Various methods which depend upon direct integration of the sensitivity equa- 
tions (ref. 9), or on finite difference methods (ref. 10) have been used in the past. 
The method of reference 11 was employed in this study. In this method, an efficient 
scheme for estimating sensitivities is employed. This scheme reduces the computa- 
tional effort required to compute the sensitivities normally associated with the 
methods used in ref. 9 and 10. 

A 

Once A0 i§ found (as the estimated value of A9^ from eq. (29)), the estimate 
of the unknown, 0, is found by incrementing 9g by A0. Because of the nature of 
this technique, the estimation process is repeated until the maximum likelihood cost 
function is minimized. 

Because of the difficulty in trying to implement model structure determination 
capability in the ML technique, it is generally necessary to know the structure of 
the model in advance. This is in general not a problem because an adequate model 
structure can be provided by the stepwise regression analysis. Though the ML tech- 
nique is computationally more difficult, it does provide asymptotically unbiased 
estimates. The ML technique is not applicable to partitioned data. 


AERODYNAMIC MODEL EQUATIONS 

Before parameter estimation and model structure determination techniques can be 
applied, a model describing the aircraft aerodynamics must be formulated. Such a 
model must not only be complete and realistic, but must also consider the estimation 
techniques used so that accurate estimates and adequate structure can be determined. 
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Assuming that the aircraft is free from any unsteady aerodynamic effects, and the 
body axis system discussed in appendix D is used, the aerodynamic forces and moments 
of the various aircraft configurations in question can be postulated as 


C = C (a) + C (a)q + C (a)6 ; (a = X, m) 

3. cL 3 a P n 

q \ 


(30) 


and as 


C = C (a) + C (a)p + C (a)r + C (a)S + C (a)6 ; (a = Y, 1, n) (31) 

a a. a a a. a a. r 

3 p r o 6 

a r 


In this formulation, the stability and control derivatives are postulated as func- 
tions of a only. C a (a) can include nonlinear functions of a. For example, C a (a) 

might comprise the terms C„ , C a a, C_ , or for a simple linear model C (a) 

d o a 

would only comprise C„ , C_ a. Furthermore, the time derivatives of a and 3 are 

a o a a 

left out of the models due to their near-linear relationships with the remaining 
variables. It is also assumed that coupling between the longitudinal and lateral 
aerodynamics is negligible. 

The actual aerodynamic models used for estimation depend on the type of data 
that is analyzed. For the small amplitude maneuvers, the aerodynamic models were 
postulated as polynomial functions in input and output variables. Large amplitude 
maneuvers were postulated as polynomial functions of the input and output variables 
along with spline functions (see ref. 12) in a. The use of data partitioning for 
data analysis was employed for use in analyzing both small and large amplitude 
maneuvers . 

The basic polynomial aerodynamic models were composed of linear stability and 
control terms along with terms which represented first and second order variations of 
each stability and control term with a. Such models were best suited for use in 
analyzing small amplitude longitudinal and lateral maneuvers. Since each small 
amplitude maneuver covers only a small range in angle of attack, the functional 
relationship between the stability and control derivatives and a can be found only 
by analyzing many maneuvers each of which covers a different range of a. 

Because of the aerodynamic nonlinearities commonly associated with large 
amplitude maneuvers, low order polynomial models may be inadequate to represent the 
aerodynamics. The nonlinearities in such aerodynamics can be approximated by 
including higher order terms in the polynomial expressions. However, such expres- 
sions often lead to large covariances on the parameter estimates and poor prediction 
properties for the model. This problem can be avoided by incorporating spline func- 
tions in the model. Spline functions (extensively discussed in ref. 12) do not have 
the problems associated with polynomial modeling because they are non-zero only for 
prescribed intervals in a. Since each spline function covers a small range of a, 
the splines defined on these intervals can approximate the nonlinearities over the 
large a range quite well using low order terms. 
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Partitioning the measured data with respect to a allows for stability and 
control parameter estimates to be found from large and small amplitude maneuvers 
without incurring the problems associated with higher order polynomial models. The 
data is partitioned with respect to a such that each partition reflects variations 
of approximately 1° or 2°. For these small variations in a, the stability and 
control parameters in each partition can be easily approximated as functions of the 
median value of a in each partition. The stability and control terms from each 
partition can then be plotted against the median values. Such plots reflect the 
functional relationships between the stability and control derivative and a. 

DISCUSSION OF RESULTS 

The flight test program on which this study was based made available statis- 
tically useful amounts of data for the determination of each specific wing geometry's 
stability and control derivatives. To be statistically useful, flight data must 
contain a large enough number of points and in addition these points must cover all 
ranges of interest in the independent variables. Such data distribution is required 
for consistent estimates of trends in and relationships between variables. Because 
of an engine failure in the test aircraft prior to the end of the planned test 
period, there were no quasi-steady flight data available for analysis for the basic 
wing geometry. Also, there were no combined maneuvers available for analysis for 
both the FOLE and the OLE geometries. Though certain types of maneuvers were not 
available, in all cases there were large amounts of both longitudinal and lateral 
data available for both stepwise regression and maximum likelihood analysis. Table 2 
lists the numbers of data points (at a sample rate of 20 points per second) available 
for analysis by maneuver type for each wing geometry. 

The nondimensional stability and control derivatives for each wing geometry were 
found using both the stepwise regression and maximum likelihood techniques. The 
majority of the available data was analyzed using the stepwise regression technique 
because the technique was significantly less expensive in terms of computer time than 
the maximum likelihood method and could be applied to partitioned data. Figure 4 
illustrates a comparison between the SR estimates and the ML estimates of the 
Cm. derivative for the SLE geometry. Tables 3 and 4 show examples of the 

typical results given by both estimation schemes. Inspection of the tabulated SR and 
ML results indicate some differences in the unknown parameter estimates. Whereas the 
SR estimates could always be found, much difficulty existed in determining the ML 
estimates wing to high correlations between the estimates. The ML estimates of weak 
parameters also tended to show large standard errors because of their low sensitivi- 
ties. 

In the analysis of the aircraft, it was assumed that the spinning propeller had 
negligible effect on the aircraft aerodynamics. However, such an assumption is not 
valid for an operating propeller which causes thrust effects. Because of such 
effects, problems can result in estimating the X-force derivatives. Hence, only the 
vertical and lateral force derivatives and the pitching, rolling and yawing moment 
derivatives were determined in this work. 

The final estimates of the stability and control derivatives for all of the wing 
geometries studied are shown in figures 3-10. These figures, which represent the 
functional relationship between the stability and control parameters and a, were 
obtained by passing a curve through the combined stepwise regression and maximum 
likelihood estimates. Where the parameters were estimated with little scatter, the 
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final stability and control parameter curves were created using numerical least 
squares curve-fitting techniques. For the parameters which had much greater scatter 
in their estimates, the final curves were obtained by hand-fitting a curve through 
the parameter estimates. 

If possible, stability and control parameters determined from flight data should 
be compared to estimates obtained from wind-tunnel tests. Since there were no tunnel 
data suitable for analysis, the final estimates of each wing geometry's stability and 
control parameters were verified by comparison to the theoretical estimates discussed 
in appendix B and by analysis of their aircraft behavior prediction capability 
(discussed in appendix C). Comparison of the theoretical estimates to the flight 
estimates and comparison between the measured and predicted time histories indicate 
that the flight data estimates were physically realistic and that the estimated aero- 
dynamic models were adequate. 

A discussion of the comparison of the stability and control parameters for the 
basic, FLE, SLE, and FOLE wing configurations follow, along with a separate compari- 
son between the similar OLE and FOLE geometries. 


Vertical Force Parameters 

Analysis of the vertical force parameters shown in figure 5 indicates signifi- 
cant differences between the SLE/FOLE and the basic/FLE wing geometries. Both the 
SLE and FOLE show greater magnitudes in their vertical force coefficients as a 
increases above 10° as compared to the basic or FLE geometries. Also, the SLE/FOLE 
configurations tend to Indicate less stalling tendency. There also exists greater 
increases in the magnitude of C for the SLE/FOLE geometries than for the basic 

or FLE geometries. Such differences may occur because both the SLE/FOLE wings 
possess irregular leading-edge geometries which generate vortices into the flow over 
the upper surface of the span. Such vortices tend to decrease the spread of separa- 
tion on the upper surface of the wing, thus reducing stalling and separated wake 
effects. 


Pitching Moment Parameters 

The pitching moment parameters are presented in figure 6. Inspection of this 
figure indicates that all the wing geometries possess static stability in pitch at a 
CG location of 31% c. Also of interest Is the fact that all the configurations 
analyzed exhibit no ability to fly at trimmed conditions without some negative 
stabilator deflection. Analysis of the other pitching-moment parameters shows 
significant differences in the pitch-damping derivative, (V , between each wing 

geometry. As the aircraft enters the stall regime, all the^geometrles studied show 

decreases in C' . However, the FOLE geometry shows a maior loss in C' as the 
m m 

stall regime is penetrated. The FOLE along with the basic and SLE geometry's pitch- 
damping become more negative as the aircraft enters the post-stall flight regime. 

The stabilator effectiveness exhibited no significant variation between any of the 
geometries studied. It is usually expected that the trends in (V should follow 

^ h 

the trends in C since the stabilator effectiveness is approximately proportional 

\ 

to the tail lift-curve slope. Since C' is also a function of the downwash slope 
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at the tail, the tail downwash must be varying to such a degree as to reduce the 

effect of the tail's lifting force on the stabilator effectiveness. In the future, 

analysis of the flow-field associated with these modifications, specifically wake- 

fuselage and wake-tail interaction, can also be of use in explaining the variation in 

trend between C 7 . and C™. 

^ 6 h m «h 


Lateral Force Parameters 

Though there were no discernable trend differences in the lateral force para- 
meters between the regular and irregular wing geometries, there were noticeable 
differences in the magnitudes of the lateral force parameters for each geometry. As 
figure 7 shows, there were variations in the Cyp derivative between each 

geometry. The basic wing showed a pronounced loss in magnitude of Cy^ , whereas 

the SLE, FLE, and FOLE configurations indicated less or no decrease in the magnitude 

of C Vo . For these four geometries, the C v - derivative generally behaved 
/p y 6 y 

as expected as the aircraft entered the stall regime, i.e., Cyg decreased in 

magnitude as angle of attack increased due to wake interaction with the tail. The 
SLE wing exhibited a drastic drop in Cyg as the aircraft left the stall 

regime. Upon full penetration into the post-stall regime, the SLE geometry then 
experienced an increase in Cyg . Though the remaining lateral force 

parameters are generally insignificant and consequently difficult to estimate, both 
the Cy r and Cyg derivatives for the SLE, FLE, and the FOLE geometries 

were found. Both the FLE and FOLE configurations had similar Cyg and Cy r 

derivatives, while the SLE geometry showed major differences in magnitude and sign of 
these derivatives. 


Rolling Moment Parameters 


The rolling moment parameters are presented in figure 8 . Though there were no 
significant variations in magnitude in these parameters between wing configuration 
type (variations less than 20 percent), two derivatives showed variations with 
respect to leading-edge regularity. The FOLE geometry showed an increase in aileron 
effectiveness in roll, Cjg , of approximately 30 percent along with an 


increased magnitude in 



in the pre-and post-stall flight regimes. 


Such 


variations in these modification's rolling characteristics are attributable to 
reduction in the separated flow over the outboard surface of each irregular 
geometry's wings. 
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Yawing Moment Parameters 


The yawing moment parameters, illustrated in figure 9, showed major variations 
with respect to wing geometry. The weathercock stability parameters, C n g , would 


in general decrease in magnitude as a increased due to wake-fuselage interaction. 
This trend is exhibited by the FOLE, FLE and basic wing configurations. The SLE 
geometry indicated increases in C n g as the aircraft approached the stall regime. 

The SLE modification then exhibited a drastic decrease in the weathercock stability 
as the post-stall regime was penetrated, whereas the FLE, FOLE, and basic geometries 
exhibited C no derivatives of consistent magnitude. Both the irregular 

leading-edge geometries possessed greater rudder effectiveness parameters, 

Cn, , than the FLE or basic wing configurations. Of special interest was the 
■o r 

basic wing which possessed total loss of rudder effectiveness as a approached 23°. 
The yaw-damping parameter, C nj _, was only identifiable in the pre-stall flight 

regime for the basic wing geometry. All the geometries studied showed loss in 
yaw damping as a increased, with total loss in yaw damping occurring for the SLE 
geometry at a = 17°, and total loss in C nr for the FLE geometry occurring at 


a = 23°. The yaw due to roll parameter * ^n n , was also not well identifiable for 

the basic wing configuration. Though the FLE and FOLE geometries showed similar 
trends in C n , the SLE configuration exhibited large variations in magnitude in 


1 p as a increased. 


Comparison of OLE and FOLE Geometries 


Figure 10 illustrates the variation in the stability and control parameters with 
a for the OLE and FOLE modified wing geometries. Because both geometries are very 
similar, major differences in the values of the stability and control derivatives for 
these geometries would not be expected, and analysis of these estimates confirms 
this. Though there were variations of approximately 15 percent between the 
parameters of the OLE and FOLE geometries, only the Cg^ derivative showed any 


major differences in magnitude (but not in trend). 


The parameter for the 


OLE geometry was 30 percent smaller than that of the FOLE configuration for the 
pre-stall regime, and there existed up to a 50 percent decrease in for the 

OLE geometry in the post-stall regime as compared to the FOLE wing. This variation 



irregular leading-edge modification’s induced vortices. 


CONCLUSIONS 

The effects of wing-leading-edge modification on a general aviation aircraft’s 
stability and control parameters have been presented. The modifications consist of 
varying the geometry of the leading-edge by addition of various drooped sections. 
Flight data were analyzed using two parameter extraction techniques for five 
different wing geometries (the basic wing and four configurations possessing 
leading-edge modifications). The analysis of the flight data resulted in complete. 
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verifiable sets of estimates for each wing geometry's stability and control 
parameters for 4° < a < 24°. From comparison of these stability and control 
parameters, the following points were noted: 

(1) The basic wing suffers from loss of lift due to stalling and the resulting 
loss of roll-damping and weathercock stability. This could result in less than 
favorable stall/spin characteristics. 

(2) The full leading-edge (FLE) modification experiences the same problems in 
lift and stability as the basic wing. The results indicate that the FLE wing 
provides no mechanism for the improvement of stall/spin characteristics. 

(3) The segmented leading-edge (SLE) wing configuration exhibits reduced 
weathercock stability and roll-damping for angles of attack greater than 20°. 

However, the SLE wing can have better stall/spin characteristics than the basic or 
FLE wings because of its improved flow over the wing surface. 

(4) The faired outboard leading-edge (FOLE) modification suffers from poor 
pitch-damping as the aircraft enters the stall regime. However, this modification 
can provide the best stall/spin characteristics of all the modifications tested 
because of its improved lift, roll-damping, and weathercock stability parameters. 

(5) The outboard leading-edge (OLE) and FOLE wing configurations possess very 
similar stability and control parameters owing to similarity in their geometries. 

In the final analysis, complete models of aircraft aerodynamics for each of the 
geometries studied were determined. A systematic study of these configurations has 
shown quantitatively that the geometries which possess irregular leading-edge 
geometries are most favorable for improving the stall/spin characteristics of general 
aviation aircraft. 
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APPENDIX A 


DATA CORRECTION EQUATIONS 

I. Upwash/sidewash correction equations for FDV sensors In degrees and knots are 


B,L 

= -1.597 + 

0.8196a u 
L , M 

(Al) 

B,R 

= -0.9003 + 

0 * 8196a R,M 

(A2) 

B,L 

= 3 l,m 


(A3) 

B,R 

= 3 r,h 


(A4) 

B,L 

= -2.757 + 

1.0449V w 
L , M 

(A5) 

B,R 

= -1.946 + 

1 * 0160v r,m 

(A6) 


II. Transformation of left and right boom velocity measurements to aircraft CG 
location 


"L = V B,L COs 3 B,L COSa B,L “ q2 L + ry L 
U R = V B,R C ° s 3 B,R C ° Sa B,R “ qZ R + ry R 
V L = V B,L Sln3 B,L " rX L + PZ L 
V R = V B,R Sln3 B,R rx R + pz R 
W L = V B,L COS0 B,L Slna B,L “ Py L + qX L 
W R = V B,R COS0 B,R Slna B,R " Py R + qX R 

y 2 2 2.1/2 

V L = (u L + V L + V 


u ( 2 . 2 . 2.1/2 

V R = (u R + V R + W R } 


(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 
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III. Transformation to obtain resultant velocity, attack and sideslip angles, and 
linear accelerations at aircraft CG location 


V = 1/2(V l + V R ) 


W T W 

a = 1/2 [tan — + tan — 1 


-1 V L . . -1 V Ri 


8=1/2 [sin -rr- + sin ] 

L R 


2 2 

a x = a^ + [(q + r )x x - (pr - r)y x - 


2 2 

ay = a Y + t(P + r )y y - (pq + r)x Y + 
M 


2 2 

a z = a z + t(q + p )z z - (pr + q)x z - 
M 


(pr - q)z x l 
(qr - p)z y ] 
(qr + p)y z ] 


(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 


IV. Calculation of the orientation of the gyro Inertial reference system based on 
accelerometer measurements for straight and level unaccelerated flight. 


<p = 0 (A21) 

u 


6 = sin *(a ) 

VI A 

U 


(A22) 




u 


sin 


-L 


(-ay /cos0 u ) 
u 


(A23) 


V. Transformation of measured attitude angles to Euler angles based on the gyro 
inertial reference system and the measured attitude angles 


0 


sin *((sin0 cos0..cos'|' w - sin<j> cos0 tI cos0 v/ sini(O 

U M M U U M M 

+ cos<{> cos0 sln0 w ) 
u u M 


(A24 ) 
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(A25) 


sin<p.,cos0 = sin0 (simp.,cos<p., - simp sin0 cosip ) 

M U M M M M M 

+ sin <p co 80 (simp.,sin0..sin<p + cosip..cosip ) 
u u M M M MM 

+ cos<p cos0 simp cos0 w 
u u u M 


cos<p w cos0 = -sin0 (cos<p 1/ cosip..sin0 v . + sinip.,sin<p 1L ,) 

M u M M M MM 

+ simp cos0 (simp ll cos<p v ,sin9., - cosip^simp.,) 

U U M M M MM 

+ COS<p COS0 COS<p,,COS0 w 
u u M M 


sin<p M cos0 

^ ~ tan simp cos0 
M 


simp w cos0 = cos<p cos0.,sinip w + simp sin0 
M u M M u M 


cosip w cos0 = cos0 cos0 w cosip + sin<ji sin0 cos0 simp 
M uMM uuMM 

- cos<p sin0 sin0 


M 


tan 


_jSimp M cos0 


cosip w cos0 

M 


(A26 ) 


(A27) 

(A28) 

(A29) 


(A30) 
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APPENDIX B 


THEORETICAL PARAMETER ESTIMATES 

Theoretical estimates of aircraft stability and control parameters were obtained 
from consideration of aircraft geometry. The knowledge of each aircraft's planform, 
wing camber, and control surface locations allowed for the determination of basic 
aerodynamic characteristics by use of vortex-lattice (VL) methodology. The vortex- 
lattice analysis was completed using a FORTRAN program (14) which determined aero- 
dynamic forces and moments by modeling the aircraft as a 3-D planform on which horse- 
shoe vortices were placed. An example of the planform representation of the aircraft 
is presented in figures ll(a-b). The longitudinal stability and control parameters 
along with one lateral stability parameter were obtained from the planform in 
figure 11(a), while the remaining lateral stability and control parameters were 
obtained from the planform represented in figure 11(b). 

Because of the nature of the vortex-lattice program, the complex planforms 
represented by each different wing-leading-edge geometry can be easily modeled. This 
was a benefit since standard empirical methods (ref. 15) cannot account for the 
irregular geometries that wing-leading-edge modifications possess. However, because 
vortex-lattice methods are based on potential flow theory, viscous effects are un- 
accounted for. Hence, lift coefficients are of larger magnitude than in actual 
flight conditions and separated wake interactions with the fuselage and tail do not 
exist. Since the theoretical estimates are only valid for attached flow (a up to 
8°), and also since the propeller's ability to augment the lift is not considered, 
the theoretical study is generally most useful in estimating the order of magnitudes 
of the various parameters. 

The basic aerodynamic characteristics computed from the vortex-lattice analysis 
were coupled with the aircraft stability and control theory presented in refer- 
ences 16-21 to arrive at a set of linear stability and control parameters. 

Owing to difficulty in modeling the proper horseshoe vortex arrangement for the 
SLE wing geometry, theoretical estimates of that geometry's longitudinal parameters 
were not obtained. Furthermore, only the lateral parameters of the basic wing con- 
figuration were found since the lateral planform cannot model the various wing- 
leading-edge geometries. 

The following formulae were used to obtain the stability and control 
parameter estimates: 


LONGITUDINAL 


The vortex-lattice program computes 
small a, then 


°L ' C m » a L=0’ °L ’ and C 1 * 
a q q p 


Hence, for 


C Z “ _C L (Bl) 

a a 



(B2) 
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(B3) 



Since the tail is a horizontal stabilator, the change in lift in the aircraft with 
respect to tail deflection equals the scaled tail lift-curve slope. Using the VL 
program to the model the tail, the tail lift-curve slope, aj is computed. Hence 



(B4) 


V "Vi < b5 > 

s h 


To account for wing-wake downwash, the downwash slope, 3e/3a, was computed from a 
theoretical method presented in reference 18. Accounting for wing-wake then results 
in the equations 



. 2v a ^ 
H a l - 3a 
c 


(B6) 




4m m* 
a 


°L 


(B7) 





+ 



(B8) 


C' 

m 6 


h 



+ 


c c 

4m Z x 

a 6 


(B9) 


LATERAL 

The lateral stability and control parameters, being based on the nonsymmetric 
lateral planform, are generally more difficult to estimate than the longitudinal 
stablity and control parameters. Because the VL program used in this study must deal 
with symmetric planform geometries, the upper and lower surfaces and vertical tail 
are each modeled separately (see fig. 11(b)). Such a method is discussed in refer- 
ence 19. The analysis of upper and lower fuselage surfaces yields 
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(BIO) 





up 


low 



2 



a 


F 


(Bll) 



s 

C = - ~ a sinr 

Y 6 S a 
a 



(B12) 


(B13) 


(B14) 


The weathercock stability parameter, C , is generally a function of Cy and the 

^0 * 8 

distance between the lateral center of pressure and the aircraft's CG location. The 
aircraft's lateral center of pressure can be approximated as the center of area of 
the lateral planforra. Completing the theoretical analysis then results in 



a a 


S_ In 

C n fi = ’ a rTT 
r 


(B15) 


(B16) 


(B17) 


(B18) 


(B19) 
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(B20) 
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(B21) 
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1 
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a 



2b<S 


a 


[C 


V 5 

0 a 





(B22) 


(B23) 


where 6 is any aileron deflection at which CL can be determined. 
a l L 0 

The theoretical estimates of each modification's stability and control deriva- 
tives are presented in tables 5 and 6. Plots illustrating the comparison between the 
flight and VL estimates of representative stability and control derivatives are pre- 
sented in figures 12(a-b). Inspection of these figures and comparison between the 
tabulated VL estimates (tables 5-6) and the flight estimates (figs. 5-10) indicate 
very good agreement in magnitude between both sets of estimates. This indicates that 
the flight data estimates are physically realistic. 
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APPENDIX C 


MODEL /PARAMETER VERIFICATION VIA PREDICTION CAPABILITY ANALYSIS 

The estimated aerodynamic models were verified by analyzing how well the models 
could predict aircraft response to known control Inputs. Knowing an adequate model, 
the prediction should match the actual aircraft response to the control inputs. This 
also is useful as a check on the estimated stability and control derivatives, since 
they are incorporated Into the aerodynamic models. 

A 6-DOF aircraft response is generally predicted by integrating equations D1-D8 
using some numerical method. This requires that the Initial condition of the states, 
the complete time history of the control Inputs, and the equations for aerodynamic 
force and moment coefficients as functions of the aircraft states and controls are 
known. This Information is directly available from the stability and control param- 
eters and from the digitized flight data. 

In order to simplify the analysis of the aircraft, It was assumed that the 
spinning propeller had negligible effect on the aircraft's aerodynamics. However, in 
actual flight conditions, the propeller is not free but spins at the idle engine 
power RPM. This causes thrust which affects the longitudinal force on the aircraft, 
which in turn affects the longitudinal velocity, u. Since the thrust effect of the 
propeller was not known, the longitudinal velocity was taken as a known state from 
the flight data. Therefore, only a 5-DOF system remained which was solved by 
integrating equations D2-D8. The integration was completed using a fourth order 
Runge-Kutta scheme. 

Examples of the predicted and measured aircraft response from this analysis are 
shown in figures 13(a-b). The measured data used in comparison to the response pre- 
dicted from the estimated models was itself not incorporated in the parameter 
estimation. As inspection of these figures indicates, the prediction capabilities of 
the estimated models are quite satisfactory, thereby verifying the adequacy of the 
estimated aerodynamic models. 
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APPENDIX D 


EQUATIONS OF MOTION 

Figures 1 and 14 illustrate the body-axis system along with 
the velocities and angular rates and the Euler angles neccessary 
motion of an aircraft. Under the following assumptions 

1. The airplane is a rigid body 

2. The effect of the spinning propeller is negligible 

3. The airplane has a plane of symmetry in the XZ plane 

4. The rotation of the earth is neglected 

the equations which explain the motion of an aircraft are: 


u = -qw + rv - g sin 0 


+ PV^S c 

2m C X 


v = - 


pV s 

•ru + pw + g cos 0 sin <j) + — — ■ 


w = 


pV s 

= -pv + qu + g cos 0 cos <j> + — ^ — 


I f ~ 1 Z l XZ, . pV 2 Sb „ 

P - qr -r— + Y-(pq - r) + -jp- C 


4 . pr III!* + Is (r 2 - P 2 > + c 

X Y X X 2I Y m 


f - P q lijlll + !“(* - , r) + ^ c 
Z X Zi z 


0 = q cos <}> - r sin 

<j) = p + (q sin <J> + r cos <j>) tan 0 


the directions of 
to describe the 


(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 
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TABLE 1.- AIRCRAFT GEOMETRY, MASS AND INERTIA CHARACTERISTICS 



Basic 

FLE 

SLE 

FOLE 

OLE 

Wing (Modified NACA 63 2 A415) 

O 

Area, m 

13.56 

13.61 

13.74 

13.59 

13.59 

Aspect ratio 

7.395 

7.318 

7.249 

7.329 

7.329 

Span, ra 

9.98 

9.98 

9.98 

9.98 

9.98 

Mean aerodynamic chord, m 

Dihedral, deg 
Root Incidence, deg 
Tip incidence, deg 
Aileron area (each), m^ 
Aileron span (each), m 
Aileron deflections, deg 

1.34 

1.39 

1.38 

+ 6.5 
+ 3. 

+ 1 • 

.94 

1.65 

-20/+10 

1.40 

1.37 

Stabilator (NACA 63AJ012) 

2 

Area, m 
Aspect ratio 
Span, m 

Mean aerodynamic chord, m 
Tail length, ra 
Tail volume 
Deflections, deg 



2.516 

4.33 

3.302 

.762 

4.383 

.607 

-15/+2 



Vertical tail (NACA 63A!012) 

2 

Area, m 
Aspect ratio 
Span, ra 

Mean aerodynamic chord, m 
Tail length, m 
Tail volume 
Rudder area, 

Rudder root chord, m 
Rudder tip chord, ra 
Rudder deflections, deg 



1.356 

1.46 

1.405 

1.023 

4.138 

.041 

.429 

.374 

.318 

+I-25 



Mass 

Flight averaged, kg 

1100 

1060 

1069 

1015 

1082 

Inertia (flight averaged) 

X X» k g~ m 2 

2245 

2351 

2406 

2319 

2380 

I Y , kg-nr 

3168 

3067 

2618 

2457 

3016 

I z , kg-nr 

4109 

4846 

4455 

4193 

4876 

ixz* k g _m2 

157 

172 

168 

180 

176 
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TABLE 2.- TABULATED 

NUMBERS 

OF DATA 

POINTS 

AVAILABLE 

FOR ANALYSTS 



Wing 

configuration 


Maneuver type 

Basic 

FLE 

SLE 

FOLE 

OLE 

Longitudinal 

5800 

7700 

4400 

6800 

7700 

Lateral 

4600 

27000 

3000 

12000 

13000 

Combined 

2600 

3300 

6300 

* 

* 

Quasi-steady 

* 

4500 

3600 

4800 

2100 


*No data available. 
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TABLE 3.- PARAMETER ESTIMATES AND THEIR STANDARD ERROR AS DETERMINED BY 
SR AND ML FOR A LONGITUDINAL MANEUVER 


FLE Wing Geometry 

(flight 185, record 13a, 0° < a < 9°) 



SR Estimates 

ML 

Estimates 

Parameter 

Value 

Standard 

error 

Value 

Cramer-Rao 
lower bound 
standard error 

C Z 

z 0 

-.22 

— 

-.62 

.002 

c z 

-3.99 

.019 

-4.23 

.079 

c z 

q 

-12.6 

.41 

-3. 

1.23 

c 7 

Zfi h 

-.62 

.022 

-.71 

.064 

C' 

m 0 

-.072 

— 

* 

* 

C' 

m a 

-.025 

.0091 

-.049 

.0042 

C' 

m 

q 

-15.1 

.20 

-16.1 

.48 

C' 

111 r 

-1.86 

.011 

-1.93 

.023 

•ff 

After second 

pass through 

program. 

C' was < 

m o 

10 - ^, therefore C' 
m i 


set to zero to decrease computational effort. 


was 
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TABLE 4.- PARAMETER ESTIMATES AND THEIR STANDARD ERROR AS DETERMINED BY 



SR AND 

ML FOR A LATERAL 

MANEUVER 



(flight 

SLE Wing Geometry 
124, record 24a, 0° < a 12 

°) 


SR Estimates 

ML 

Estimates 

Parameter 

Value 

Standard 

error 

Value 

Cramer-Rao 
lower bound 
standard error 

Cy 

Y e 

-.34 

.005 

-.35 

.008 

Cy 

r 

-.80 

.044 

* 

* 


.17 

.012 

-.04 

.017 


-.062 

.0077 

-.04 

.016 

% 

-.084 

.0015 

-.085 

.0056 

% 

-.390 

.0067 

-.494 

.0033 


.32 

.013 

.14 

.004 

C, 

X «a 

-.182 

.0035 

-.19* ** 

— 

Ci 

.024 

.0022 

.016 

.0008 


*Poor or non-existent estimates due to low sensitivities. 

**Due to high correlation with other parameters, a fixed estimate was 
supplied based on SR estimates. 
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TABLE 4.- Concluded 



SR 

Estimates 

ML 

Estimates 

Parameter 

Value 

Standard 

error 

Value 

Cramer-Rao 
lower bound 
standard error 

% 

.020 

.0016 

.035 

.0003 

°”p 

-.157 

.0070 

-.056 

.0024 

C„ 

n r 

-.14 

.013 

-.139 

.0013 

c % 

* 

* 

.0023 

.00075 

H 

-.047 

.0023 

-.046 

.0004 

*No estimate 

obtained 

due to Insignificance of 

this parameter. 


32 



TABLE 5.- THEORETICAL ESTIMATES OF LONGITUDINAL PARAMETERS (31% c) 


Parameter 

Basic 

FLE 

OLE 

c 7 

Z 0 

-.444 

-.434 

-.508 

iX 

-5.394 

-5.474 

-5.403 

CZ q 

-11.427 

-11.472 

-11.282 


-.697 

-.697 

-.697 

C' 

m 

a 

-.304 

-.118 

-.259 

C' 

m 

q 

-22.237 

-22.170 

-22.132 

C' 

-2.288 

-2.288 

-2.288 
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TABLE 6.- THEORETICAL LATERAL PARAMETER ESTIMATES FOR BASIC 

WING CONFIGURATIONS 


Parameter Value 

Parameter 

° y b 

.477 

% 

% 

* 

\ 

% 

.396 

c i 

r 

^a 

.034 

X 

° y 6 

r 

.143 

c l 

r 

♦Generally 

Insignificant. 


Value 

Parameter 

Value 

-.096 

% 

.070 

-.540 

s 

-.065 

.079 

% 

-.164 

-.183 

X 

.037 

.029 


-.059 
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Small amplitude SR estimate 
Large amplitude SR estimate 
Partitioned SR estimate 
Spline SR estimate 
ML estimate 


A 
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7.- Lateral Force Parameter Estimates 




Figure 7.- Concluded 
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